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I/") ■ Abstract 

(N ; 

Context. In the framework of the study of supernova remnants and their complex interaction with the interstellar medium, we report on an 
XMM-Newton EPIC observation of the Galactic supernova remnant W49B. 
' Aims. We investigate the spatial distribution of the chemical and physical properties of the plasma, so as to get important constraints on the 
, physical scenario, on the dynamics of the supernova explosion, and on the interaction of the supernova remnant with the ambient interstellar 
O |' clouds. 

' Methods. We present line images, equivalent width maps, and a spatially resolved spectral analysis of a set of homogeneous regions. 
»j . Results. The X-ray spectrum of W49B is characterized by strong K emission lines from Si, S, Ar, Ca and Fe. In all the regions studied, the X-ray 
^ spectrum is dominated by the emission from the ejecta and there is no indication of radial stratification of the elements. A high overabundance 
of Ni (Ni/Ni = 10+j) is required in the bright central region and the previous detection of Cr and Mn line emission is confirmed. Spectra 
• • ■ are well described by two thermal components in collisional ionization equilibrium. We observe spatial variations in the temperature, with the 
. ' highest temperature found in the east and the lowest in the west. 

, Conclusions. Our results support a scenario in which the X-ray emission comes from ejecta interacting with a dense belt of ambient material, 
■ but another possibility is that the remnant is the result of an asymmetric bipolar explosion with the eastern jet being hotter and more Fe-rich 
C3 ' than the western jet. The eastern jet is confined by interaction with ambient molecular clouds. Comparison of the observed abundances with 
yields for hypernova and supernova nucleosynthesis does not directly support the association of W49B with a y-ray burst, although it remains 
possible. 

Key words. X-rays: ISM - ISM: supernova remnants - ISM: individual object: W49B 



1. Introduction 

Supernova remnants (SNRs) are powerful sources of heavy ele- 
ments and energy in the interstellar medium (ISM). The super- 
nova supersonic material (ejecta) is expelled into the ambient 
medium and produces a shock front that heats and ionizes the 
ISM. As this shock expands into the ISM, a reverse shock asso- 
ciated with the deceleration of the ejecta propagates inward (in 
the reference frame of the expanding ejecta), thus heating the 
expelled material to X-ray emitting temperatures. The resulting 
X-ray emission provides important information about the emit- 
ting conditions and the chemical composition of the ejecta. 

The galactic supernova remnant W49B is one of the bright- 
est ejecta-dominated SNRs observed in X-rays. It has a center- 
Send offprint requests to: M. Miceli, 
e-mail: miceli@astropa . unipa . it 



filled morphology in X-r ays and a shel l-like morphology in the 
radio band. According to Rh o & Petrel i 1998b . W49B is similar 
to mixed-morphology supernova remnants, but its classifica- 
tion is uncertain. W49B is one of the 10 g alactic sources with 
the h ighest surface brightness at 1 GHz jMoffett & Reynolds 
1994). The diameter of the radio shell is ~ 4' and the dis- 

tance to the remnant, estimated on the basis of HI absorption 
i (I — — ~j 

( Radhakrishnan et al. 1972), is ~ 8 kpc (considering the cor- 
rections bv lMof fett & Revnol dsl 19941) . 

In the X-ray ban d, W49B was first detected with the 
Einstein Observatory dPve etalJ ll984). Th e discovery of in- 
tense Fe XXV K line emission by EXOSAT dSrnifh et al J 19851) 
indicated that the X-ray emission of W49B is thermal and 
dominated by the ejecta. This result h as been further con- 
firmed by analysis o f the ASCA spectrum ( Fuiimo to et al 119951 
iHwang et alJ l2000). Ejecta-dominated remnants are expected 
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to be relatively young, and while the age of W49B has not yet 
been determined accurately, published age estimates are in the 
range of ~ 1000 - 4000 yr llPve et all 1 9841 ISmith et all 19851 
iHwang etalJ EoOO). The images of Si, S and Fe K emission 
lines obtained with ASCA indicate that the iron emission is 
confined to the inner part of the remnant, while silicon and 
sulfur emission appears to originate from the outer regions. 
This result has been interpr eted as evidence for a stratified 
distribution of the elements (Fui imoto et alJll995l). The anal- 
ysis of the same ASCA data carried out bv lHwang et alJ (|2000) 
showed that broadband modeling of the remnant's global spec- 
trum requires two thermal components and significant over- 
abundances of Si, S, Ar, Ca, Fe with respect to the solar val- 
ues. Recently, the comparison between data obtained in X- 
rays with the Chandra satellite and in the infrared with the 
Palomar obser vatory have revealed new details of the morphol- 
ogy of W49B iKeohane et alJ2004l . In particular, the presence 
of an elongated structure in the Fe K emission (interpreted by 
the authors as a jet), and the nondetection of a neutron star 
associated with the remnant, have been considered to make 
W49 B consistent with t heoretical predictions for a collapsar 
(e.g JMacFadven & Woosievlll999h and thus a candidate for a 
gamma-ray burst remnant. 

Here we present the analysis of an XMM-Newton EPIC ob- 
servation of W49B. The high sensitivity of XMM-Newton pro- 
vides deep insights into the spatial distribution of the elements 
and the physical conditions in the remnant which lead to im- 
portant constraints on the physical scenario. The paper is orga- 
nized as follows: Sect. |2]presents the EPIC data and describes 
the data processing procedure; Sect. [3] presents the X-ray re- 
sults in terms of the global spectrum (13. 1> . imaging (13. 2> and 
spatially resolved spectral analysis (13.31 . Sect. |4] discusses the 
results in the context of observations at other wavelengths and 
SN models. Finally, we draw our conclusions in Sect.|5] 

2. Data processing 

The XMM-Newton data presented here consist of two ob- 
servations Observation ID 0084100401 and 0084100501 (PI 
A. Decourchelle) performe d on 3 April 20 04 and 5 April 
2004 with the EPIC MOS dTurner et alJ 1200 ill and EPIC pn 
JStriider et all2 001 ) cameras, as summarized in Table^. Both 
observations have pointing coordinates a (2000) = 19 A ll m 13 s 
and 6 (2000) = 9°7'1". A third data set 0084100601 (with the 
same pointing coordinates) is heavily contaminated by proton 
flares and was not used for our analysis. 

The data were processed using the Science Analysis 
System (SAS V6.1.0). To create light curves, images and spec- 
tra, we selected events with PATTERNS 12 for the MOS cam- 
eras, PATTERNS 4 for the pn camera, and FLAG= for 
both. To eliminate the contamination by soft proton flares, we 
screened the data by applying a count-rate limit on the light 
curves at high energies (10 - 12 keV for MOS and 12 - 14 
keV for the pn camera). With a 5t — 100 s binning for the light 
curves, we set the count-rate limit at 18 counts per bin for the 
MOS cameras, and 22 counts per bin for the pn camera. The ex- 
posure times are given in Tabled For each camera we merged 



Table 1. Relevant information about the data. 



OBSJD 


CAMERA 


tap (ks)* 


Mode 


Filter 


0084100401 


MOS1 


18.6/15.8 


Large Wind. 


medium 


0084100401 


MOS2 


18.7/16.0 


Large Wind. 


medium 


0084100401 


pn 


17.0/12.9 


Full Frame 


medium 


0084100501 


MOS1 


18.7/16.8 


Large Wind. 


medium 


0084100501 


MOS2 


18.7/16.3 


Large Wind. 


medium 


0084100501 


pn 


17.0/14.6 


Full Frame 


medium 


* Unscreened/Screened exposure time. 



the screened event files of the two observations using the task 
MERGE. 

All the images presented here are superpositions of the 
MOS1, MOS2 and pn images (obtained using the E MOSAIC 
task) and are background-subtracted, vignetting corrected and 
adaptively smoothed. Despite being relatively small in ex- 
tent compared to the XMM-Newton field of view, W49B has 
a high X-ray brightness and contaminates the local back- 
ground. For this reason, we used the high signal-to-noise back- 
ground event files El_fm0000JvIl.fits, El_00fm00_M2.fits, 
El_0000fm_PN.fit s (obtained using the medium filter) de- 
scribed in detail in lRead & PonmanM2003l) . The exposure and 
vignetting correction were performed by dividing the super- 
posed count images by the corresponding superposed exposure 
maps (produced using the task EEXPMAP). The pn exposure 
maps were scaled by the ratio of the pn/MOS effective areas 
to make MOS -equivalent superposed count rate images. We 
smoothed the images adaptively to a signal-to-noise ratio of 10 
using the task AS MOOTH. The background-subtracted count 
images and the background- and continuum-subtracted line im- 
ages (Sect. 13.21 were smoothed while accounting for their vari- 
ances. The resulting template was then applied to smooth the 
associated exposure maps before performing the image divi- 
sion. The continuum under the lines was estimated by scaling 
a continuum band adjacent to the line emission using the con- 
tinuum spectral slope derived from the global remnant spec- 
trum. The energy ranges chosen for the lines and the continua 
are listed in Table |2] In all the line bands, the background is 
^ 0.5% of the total signal, except for the 7.65 - 8.0 keV band 
where it is ~ 3.6% of the total signal. In the continuum band 
4.4 - 6.2 keV, the background amounts to ~ 1.3%. In all the 
images North is up and East is on the left. 

For spectral analysis, we processed the e vent files with 
the EVIGW EIGHT task jArnaud et alJl200lh to correct for 
vignetting effects. For the background subtraction we used 
the same background event files used for the images. Spectra 
were rebinned to achieve a signal-to-noise ratio per bin > 
5cr, and the spectral fits were performed simultaneously on 
the MOS1, M OS2 and pn s pectra in the 1 - 9 keV band 
using XSPEC (Arnaud 1996). The chosen on-axis response 
files are ml,2_medv9q20t5r6 _all_15.rsp for MOS 1,2 spectra 
and epnJF20_sdY9_medium.rsp for pn spectra. All the re- 
p orted errors are at 90% confidence, as calculated according 
to lLampton et a i1 dl976l) . 
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3. X-ray emission 

The broadband (1-9 keV) X-ray morphology of W49B ob- 
tained with the EPIC cameras is shown in Fig. ^ (left panel). 
For c omparison, the VLA radio image obtained at 327 MHz 
bv lLacev et al.l (1200 II) is also presented in Fig. ^ (right panel). 
The X-ray surface brightness of W49B has its maximum at the 
center of the remnant, and shows an elongated barrel-shaped 
structure extending for ~ 5 pc along its east-west axis (labelled 
b, see Fig.Q . 

At the eastern end of the remnant there is a bright elongated 
X-ray structure (labelled a, see Fig.Q, running perpendicular 
to the X-ray barrel and parallel to the radio shell (labelled a, 
see Fig. The eastern side of W49B is bounded by a large 
molecular cloud detected in the 2.12 /urn molecular hydrogen 
narrow band image of W49B (Fig.|2l. The X-ray barrel b lies 
between two radio structures (J3 and y) which are parallel to 
the radio shell (a) and to the H2 emission. The positions of 
these structures (J3 and y) are almost coincident with two coax- 
ial ring-like structures observed in the 1 .64 yum [Fe II] image of 
W49B (Fig. The X-ray structures a (east) and c (west) lie 
outside these radio structures. No intense X-ray emission has 
been detected in the south-western region of the shell, where 
copious radio and [ Fe III emissi on are observed (5, see Fig.Q. 
According to lLacev et al . (2001), there is evidence for stronger 
interstellar absorption in this direction. It is reasonable to con- 
clude, however, that no hard X-ray emission is present in this 
region. 

3.1. Global spectrum 

The impressive global EPIC spectrum of W49B in the 1 -9 keV 
energy band is shown in Fig. [3] The spectrum is dominated by 
intense emission lines of He-like and H-like ions of Si, S, Ar 
and Ca and by a huge Fe K emission line blend. A phenomeno- 
logical modelling of the global spectrum (2 bremsstrahlung for 
the continuum + narrow gaussians for the lines) allows us to 
derive the line energies reported in Fig. [5] The emission fea- 
ture observed at ~ 7.8 keV corresponds to a possible blending 
of Fe XXV and Ni XXVII K emission lines and is discussed 
in sect. 13.41 The continuum above 2.3 keV is well described 
by a thermal bremsstrahlung model with kT = 1 .69 ke V and 
Nh = 5.3 x 10 22 cm~ 2 ; this model has been used to esti- 
mate the continuum under the lines given in Table|2] Unlike the 
AS CA spectrum fawang et all2000l) . the XMM-Newton global 
spectrum cannot be correctly fitted by a two-component ther- 
mal model (reduced^ 2 ~ 2.4). This indicates that the emission 
originates from physically inhomogeneous regions. Spatially 
resolved spectral analysis of homogeneous regions is required 
to correctly reproduce the spectra with 2-component models, 
as discussed in Sect. 13.31 

In the analysis of the ASCA global spectrum of W49B, 
iHwang et alJ J2000l) found evidence for the presence of two 
new emission line features (at ~ 5.7 keV and ~ 6.2 keV) 
that are not included in the currently available atomic database 
for X-ray spectral modeling. These were attributed to the He- 
like Cr transitions (forbidden line at 5.655 keV; resonance 
line at 5.682 keV) and the He-like Mn transitions (forbidden 



fine at 6.151 keV; resonance line at 6.181 keV). Our XMM- 
Newton data confirm this detection. Modelling the global spec- 
trum of the remnant in the 4.4 - 6.4 keV band with a ther- 
mal bremsstrahlung continuum, we found a clear excess in 
the residuals corresponding to these transitions. These lines 
can be reproduced by adding two gaussian components to the 
bremsstrahlung model to determine the line energy and flux. 
We found, respectively for Cr and Mn, line centroid energy 
E Cr = 5.66 + 0.01 keV and flux F Cr = 2.5 ± 0.4 x 10~ 5 pho- 
tons cnr 2 s _1 , E M „ = 6.19+Jj-^ keV and F M „ = 1.0 + 0.3 x 
10" 5 photons c m" 2 s" 1 , in quite goo d agreement with the val- 
ues reported bv lHwang etall fcOOOh : E Cr = 5.69+™ 2 keV and 
F Cr = 3.0^1 x 10~ 5 photons cm" 2 s _1 ; E Mn = 6.17 ± 0.05 keV 
and Fmh = 1.3^i*g x 10~ 5 photons cirT 2 s _1 . The addition of 
the Cr gaussian leads to a reduction of x 2 , fy 1 = 111 for 616 
d. o. f.; adding the Mn gaussian gives an additional 5'x 1 = 28 
(to be compared with the values obtained with the AS CA spec- 
trum: 6x 2 = 20 for 61 d. o. f. and 5'x 1 = 10). We found no 
indication for the presence of the Ti Hea line (at 4.97 keV), 
obtaining an upper limit for the flux Fj, < 3.6 x 10~ 6 pho- 
tons cirT 2 s , that is more stringent that th e one reported from 
ASCA (F Ti < 2 x 10~ 5 photons crrT 2 s^. lHwang etalfeO OO). 
The flux of the Cr and Mn lines contributes ~ 3% of the total 
flux in the 4.4 - 6.4 band. 

3.2. Imaging 

The spatially resolved spectral study of the emission in the dif- 
ferent lines provides information about the spatial distribution 
of the synthesized elements and the chemical composition of 
the X-ray emitting plasma. On t he basis of low angu lar reso- 
lution ASCA images of W49B, Fujimoto et al ( 1995) argued 
that the sulfur emission, unlike the iron emission, was associ- 
ated with the outer envelope of the remnant and not related to 
the bright 4-6 keV continuum emission from the center. 

To investigate both this claim and the suggested radial 
stratification of the elements, we produced background- and 
continuum-subtracted line images of S, Ar, Ca, Fe and Ni lines 
(see Table|5Jl. Si emission line images, shown in Fig.|6j are dis- 
cussed in Sect. 13. 3. 11 The S H-like line image is very similar to 
the S He-like line image and therefore is not shown in the fig- 
ure. All the line images present a bright centrally barrel-shaped 
morphology similar to the 4.4 - 6.2 keV continuum emission, 
as shown in Fig.0] Overall, there are no significant differences 
between the morphologies in the Si, S, Ar, Ca and Fe lines (al- 
though some detailed structures may differ as we will discuss 
later). This result is not in agreement with the line imaging re- 
sults obtained for the lower- Z elements with AS CA data (see 
Fig. 4 in Fuiim oto et al.ll 995). 

A detailed analysis of the Fe K line image in Fig.|4](lower 
left panel) shows that the iron emission differs from that lower- 
Z elements in the western part of the remnant. In particular, 
regions 6 and 7 do not have significant Fe line emission, while 
they do have significant emission from Si, S, Ar and Ca. The 
lower statistic Fe XXV+Ni XXVII image of Fig.0](lower right 
panel) presents a morphology that is fully consistent with that 
of the Fe K image. 
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Table 2. Parameters used for the production of the line images of Fig.0]and the equivalent width maps of Fig. [5] 



Emission line Line energy band (keV) Continuum energy band (keV) 

SXVKa 2.35 -2.53 2.75 - 3.0 

S XVI Ka 2.58 - 2.7 2.75 - 3.0 

Ar XVII Ka 3.05 - 3.2 2.75 - 3.0 

CaXIXKa 3.8 -4.0 3.45 -3.7 

Fe XXV Ka 6.45 - 6.9 4.4 - 6.2 

Fe XXV+Ni XXVII K 7.65 - 8.0 7.05 - 7.5 



To disentangle higher element abundances in the bright line 
image structures from higher emission measure (which might 
be caused by either higher density or higher emission volume), 
we produced equivalent width maps as shown in Fig. [5] These 
depend linearly on the abundances, but are also affected by 
the temperature, column density and ionization age. These im- 
ages have been constructed by dividing the background- and 
continuum-subtracted line images by the corresponding under- 
lying continuum. Before the division, the estimated continuum 
map was adaptively smoothed using the same template used for 
the background- and continuum-subtracted line images. As ex- 
plained in Sect. |2 the underlying continuum is estimated from 
continuum regions adjacent to the line emission, and assumes 
that the continuum has a uniform spectral slope over the en- 
tire SNR. As shown in Fig. [5] the S map is clearly different 
from the other maps, but the equivalent width images of Ar, 
Ca, Fe and Fe+Ni all show similar features. In particular, all 
these maps show higher equivalent widths in the central and 
eastern parts of the remnant than in the western part, with max- 
ima corresponding approximately to regions 1 and 3. This is 
particularly true for iron. In the Ar map there is an indication 
of non-negligible equivalent widths in the western regions 6 
and 7. 

As the S H-like and He-like line images are very similar, we 
used the combined 2.35-2.7 keV band to improve the statistics 
of the sulfur equivalent width map as shown in Fig.|5](upper left 
panel). The map appears patchy and shows no clear structures. 
Given the relatively low line energies, the map may be affected 
by local variations of Nh, as we shall see in Sect. 13. 3. T1 

In conclusion, the line emission maps of all the observed 
elements show a centrally elongated bright structure in good 
agreement with the morphology of the continuum. The equiv- 
alent width images indicate that Ar, Ca and Fe may be more 
abundant in the central and eastern parts of the remnant than 
in the west. The high values of the equivalent width are there- 
fore not restricted to the barrel-like structure, but extend to the 
eastern side of the remnant, as illustrated in Fig. |5] (lower right 
panel), where we show the contours of the iron line image su- 
perimposed on the iron equivalent width map. A clear drop-off 
of the equivalent width is observed towards the western side. 
This effect is particularly strong for iron, and less pronounced 
for the lower-Z elements which show significant equivalent 
widths in the west. While we thus find an East/West abundance 
anisotropy, we find no clear evidence for radial stratification of 
the elements. 



3.3. Spatially resolved spectral analysis 

To quantify the equivalent width map results for the element 
distributions, we performed a spatially resolved spectral analy- 
sis that allowed us to study in detail the physical and chemical 
characteristics of the plasma in different parts of the remnant. In 
particular, we analyzed the spectra extracted from seven fairly 
homogeneous regions shown in Fig. 0] (lower left panel) and in 
Fig. [5] While this set of regions does not cover the entire rem- 
nant, it does allow us to study the spatial variations of the spec- 
tral properties in the X-ray emitting plasma. In order to ana- 
lyze spectra extracted from physically uniform regions, we pro- 
duced a MOS mean photon energy map (i.e., an image where 
each pixel holds the mean energy of the detected MOS pho- 
tons) in the continuum band 4.4 - 6.2 keV. The shape and size 
of the seven regions were chosen in order to have very small 
fluctuations of the mean photon energy (~ 1%) and a sufficient 
number of photons. 

Regions 1 and 2 are located on the X-ray eastern elongated 
structure along the radio shell (labelled a/a in Fig.^, respec- 
tively at its northern and southern ends. They both exhibit high 
values in the equivalent width maps. 

Regions 3 and 4 cover the brightest part of the centrally 
barrel-shaped X-ray structure: region 3 corresponds to the max- 
imum of the continuum emission in the 4.4-6.2 keV band (Fig. 
0] upper left panel), region 4 corresponds to the highest Fe K 
line emission region (Fig. [4] lower left panel). 

Region 5 follows the western radio structure y (Fig. 
Region 6 covers the western X-ray structure c and, finally, re- 
gion 7 lies on the faint south-western structure that is seen in 
the continuum and K line emission of S, Ar, Ca, but not of iron. 

3.3.1. Spectral modelling and results 

We first modelled the extracted spectra using XSPEC 
with a single MEKAL model of an opti cally-thin thermal 
plasma in collisional ionization equilib rium ( M ewe et alJ l985. 
iMewe et alJl986llLiedahl et alll995l) . Good fits to the spectra 
require non-solar abundances for Si, S, Ar, Ca and Fe. All other 
element abundances were fixed at their solar values. To account 
for the contribution from the Cr line (Sect. l3.lt . we added a nar- 
row gaussian component with fixed energy; it was not neces- 
sary to account for the Mn line. We also added a systematic 5% 
error term to reflect the estimated 5-10% uncertainties in the 
calibration of the instrumental effective area. The reduced x 1 
values obtained using the two MOS and the pn cameras range 
from 1.3 (with 354 d.o.f.) for region 5, to 1.7 (445 d.o.f.) for 
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region 2. Above 2 keV, the models describe the continuum and 
emission line spectra quite well, but these single-temperature 
models severely underpredict the Si XIII Ka line, as illustrated 
in Fig. (left panel). The best-fit parameters for these fits are 
thus not reported he re. The same problem was encountered by 
iHwang et alJ f2000) in the modelling of the AS CA global spec- 
trum. We found that the quality of the fits does not significantly 
improve if a single non-equilibrium ionization (NEI) model is 
used instead. 

We obtained a good description of the spectra by adding 
a second collisional ionization equilibrium thermal MEKAL 
component. We require the two components to have the same 
abundances and to be absorbed by the same interstellar column 
density. The spectra are also well-described if the soft com- 
ponent has solar abundances and is in nonequilibrium ioniza- 
tion (PSHOCK model). In this case the temperatures of the two 
components (NEI+CIE) are comparable to those obtained for 
the two CIE components, but the fits give somewhat higher chi- 
squared values. For example, in region 2: Xnei+cie ~ 529.7 
(with 442 d. o. f.) and x\ C ie = 506 - 7 ( 443 ) ; in re S ion 4: 
Xnei+cie = 821 - 2 ( 673 ) and xlciE = 793 - 5 ( 674 ); in re g ion 

6 -X 2 NE i + cie = i 064 - 1 ( 843 ) ™&xla E = 10593 (W We 
therefore adopt the two CIE component model hereafter, while 

cautioning the reader that the high interstellar absorption com- 
plicates the description of the soft component and hence the 
discrimination between these two models. 

Fig.0shows a representative spectrum, extracted from re- 
gion 2, with best-fit model and residuals for two thermal com- 
ponents (right panel), and for comparison, a single thermal 
component (left panel). The best-fit parameters for all seven 
spectral regions are given in Table [3] 

The values of A^h in the different spectral regions range ap- 
proximately from 4.2 x 10 22 cm" 2 to 4.9 X 10 22 cm" 2 (com- 
pa red with the globa l Nn = 5.0^| x 10 22 cm" 2 estimated 
bv lHwang et aljEoOOh . To produce the line and the equivalent 
width images (see Sect. 13. 21 . we estimate the continuum con- 
tribution under the lines by assuming a uniform absorbed slope 
of the continuum based on the global spectrum of W49B. For 
low energy lines, spatial variation of the column density over 
the remnant can thus introduce errors in the line and equivalent 
width images. In particular, where A^h is higher than its global 
averaged value, we may overestimate the continuum under the 
line and hence underestimate the equivalent width. We observe 
this kind of effect in the S equivalent width image (Fig. [5] upper 
left panel), where there is a minimum in region 3 with a mean 
value of 0.32 keV, corresponding to the highest value of the 
column (A^h = 4.8*2 i), and a higher mean equivalent width of 
0.4 keV in region 2 which has a lower column (A^h = 4.22^'^). 

Variations in column density may also affect the Si images, 
but in this case there is no a clear relation between such vari- 
ations and inhomogeneities in the equivalent width map. This 
probably indicates that other effects are also present. As ex- 
plained above, a single temperature model provides a good 
overall description of the spectra in the 1-9 keV band, but 
fails to reproduce the observed Si line flux in the narrow energy 
range 1.7 - 1.9 ke V. It is therefore possible that the Si He-like 
emission originates predominantly in a plasma with different 
physical conditions than that responsible for the other emis- 



sion lines. As shown in Fig. [6] the He-like line morphology 
resembles that of the Si H-like emission as well as the S, Ar 
and Ca line emission. The Si equivalent width maps (not pre- 
sented here) have irregular patchy features at odds both with 
those seen for the other elements, and with the spatially re- 
solved spectral analysis. This is probably because the contin- 
uum of the hard component is not negligible at low energies, 
even though the Si He-like line emission itself is associated 
with the soft thermal component (see the right panel of Fig. Q. 

We observe larger variations for the temperature of the hot- 
ter component kTi (1.75 - 3.3 keV) than for the cooler compo- 
nent kT\ (0.7 - 1 .05 keV), but do not measure any radial tem- 
perature structure (see Tabled. The temperature of the hotter 
component is at a maximum in Region 1 at kTi, and decreases 
significantly moving southward along the X-ray elongated a 
structure to region 2. This temperature also decreases towards 
the western side of the remnant (from region 3 to region 6). 
We observe that the difference between the continuum surface 
brightness of regions 3 (i. e. the brightest region in the contin- 
uum band 4.4 - 6.2 keV) and region 4 is caused by the higher 
temperature in region 3, and not by differences in the emission 
measure. 

In Table |5J the parameter EM indicates the values of 
the emission measure per unit area, which is the product of 
the square of the particle density times the extension of the 
plasma along the line of sight. These values are in the range 
0.3 - 1.7 x 10 20 cm" 5 for the cooler component, and 0.4 - 
3.5 x 10 20 cm" 5 for the hotter component. In region 6 we have 
the highest values for both EM\ and EM^. Associated with the 
decreasing trend of Ti, we note that in region 6 the correspond- 
ing emission measure EM2 is slightly higher than in regions 3 
and 4. 

We used the best-fit values of temperature and emission 
measure to derive the plasma density (indicated with n in Table 
|5l a nd the filling facto r (f) of each component, obtained as in 
Boc chino et alJ dl~999h . These quantities were derived assum- 
ing pressure equilibrium between the two components and as- 
suming that in each spectral region the extension of the X-ray 
emitting plasma is equal to the length of the chord intercepted 
by the radio shell (approximated as a sphere with angular radius 
<p = 2.15') along the line of sight. The filling factor of the cool 
component is ^ 0.1, and the densities of the two components 
are almost uniform across the remnant. 

We found significant element overabundances in all the 
spectral regions: all the elements have fitted abundances > 1.5 
(at 90% confidence level) with respect to their solar value. This 
attests to the origin of the emission in the ejecta. Moreover, 
higher abundance values are found for the central and east- 
ern regions compared to the western regions, in good agree- 
ment with the results of our image analysis. All abundances 
are lower in regions 5, 6, and 7, with the largest decrease ob- 
served for the Fe abundance, as expected. For example, the ra- 
tio (Ar/Ar Q )/(Fe/Fe G ) is 0.7 ± 0.1 in region 4 and 1.8 ± 0.3 
in region 6, and the ratio (Ca / 'Ca ) / '(F 'e / F e ) is 1.2 + 0.2 in 
region 4 and 1.8*q| in region 6. This is in agreement with the 
considerations discussed at the end of Sect. 13.21 where regions 
6 and 7 show low equivalent widths for iron, but nonnegligi- 
ble line emission and line equivalent widths for the lower-Z 
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elements. All the abundances reported in Tableware given as 
the ratio of the number density of the respective element to 
the number density of hydrogen expressed in units of the so- 
lar value of this ratio. In a pure-ejecta plasma (where hydrogen 
emission is negligible) this ratio would be ill-defined, and the 
spectral shape would depend only on the relative abundances of 
the heavy elements. We verified that this is not the case for our 
fits. In fact if we scale all the abundances by the same factor / 
we obtain significantly higher^- 2 values for / ^ 1.5. Therefore 
the abundances in Tableware well defined and above solar. 

In each spectral region we also obtained another x 1 mini- 
mum with kT\ ~ 0.2 keV. This value is in agre ement with the 
results obtaine d with the A S CA data by Hwa ng et alJ d2000l) 
and iKawasaki et alJ J2005h . In this case, however, we found 
very high values for the emission measure per unit area of the 
cool component (EM\ ^ 10 22 cm 5 ), which would imply un- 
realistically high values for the filling factor (f\ ~ 0.5) and 
the density (n\ ~ 50 cm -3 ). In what follows, we will therefore 
discuss only the results presented in Table|5] 

3.4. Evidence for Ni overabundance 

Thanks to the high XMM-Newton effective area at high energies 
it was possible to detect a strong emission feature at ~ 7.8 keV, 
with possible contributions from both Fe XXV and Ni XXVII. 
The Ni XXVII emission line has been previously detected in 
only i n one SNR, Cas A dBleeker et al.ll200lllWillingale et alJ 
2002). This emission blend can be used to set constraints on 
the Ni abundance in W49B. Since the photon statistics for the 
seven individual spectra do not allow us to resolve the features 
in the 7.65 - 8.0 keV band in detail, we extracted the spectrum 
from the larger elliptical region shown in Fig. |4] (lower right 
panel). This region covers the area of highest surface bright- 
ness in the 7.65 - 8.0 keV band. We fitted the spectrum in the 
1-9 keV band with the 2-component thermal model used for 
the spatially resolved spectral analysis. Table|4]gives the best- 
fit parameters derived both by fixing the Ni abundance to the 
solar value (second column) and by fitting the Ni abundance as 
a free parameter (third column). We can better reproduce the Fe 
XXV and Ni XXVII K line blends by fitting for the Ni abun- 
dance. The F-test gives the probability that the improvement 
of the fit is insignificant as < 10~ 10 . We conclude that there is 
convincing evidence for Ni overabundance (Ni/Nio = 10^ 2 ) in 
the centrally barrel-shaped structure of W49B. 



3.5. Overionization 

Through their analy sis of the ASCA global spectrum, 
IKawasaki et al.l (|2£>05) claimed the presence of "overionized" 
plasma in W49B. In particular, they measured the intensity ra- 
tio of the H-like to the He-like Ka lines of Ar and Ca to ob- 
tain the ionization temperature T z , and found that T z is higher 
than the electron temperature T e estimated from the continuum 
of the X-ray spectrum. We performed the same analysis by 
modelling the XMM-Newton global spectrum in the 2.75 - 6.0 
keV band. We model the continuum and estimate T e with a 
zero metallicity VMEKAL component and include four narrow 



Table 4. Best-fit parameters for the spectrum extracted from 
the elliptical region shown in Fig.@](lower right panel) and re- 
produced by two thermal components in collisional ionization 
equilibrium. All errors are at the 90% confidence level. 



Parameter 


Ni/Ni Q = 1 


Ni thawed 


jV H (10 22 cm- 2 ) 


4.60*}* 


4.68!°- 


kT { (keV) 


0.87!g;« 


0.863£ 


em; (10 20 cm" 5 ) 


i i q+0.04 
1-1 -0.02 




kT 2 (keV) 


9 tc+0.02 
J -0.03 


9 9+0.2 
■^-O.l 


EM\ (10 20 cm- 5 ) 


i 7Q+0.06 
1 • ' 7 -0.04 


1 09+O.O6 


Si/Si 


T "2+0.1 

J - -0.3 


3.3 ± 0.2 


s/s e 


3.7 ±0.1 


Q 7+0.1 
•'• -0.2 


Ar/Ar e 


4 9+0-4 


4 7+0.3 


Ca/Ca s 


6.3 ±0.5 


6.4 ± 0.5 


Fe/Fe 


5 5 +01 


6 +al 

"•"-0.2 


Ni/Ni Q 


1 


10^ 


x 2 /d. o. f. 


1847.8/1327 


1584.3/1326 



Emission measure per unit area. 



gaussian components to model the Ar and Ca lines and a fifth 
gaussian component for the Cr line. We considered absorption 
column densities of N H = 4 6x 10 22 cm" 2 and N H = 5.23 xlO 22 
cirr 2 (the value adopted by Kawasaki et alJ J2005I) '). For those 
column densities, we obtained electron temperatures of kT e = 
1 .80 + 0.03 keV and 1 .66!°;° 2 keV, respectively. We find no ev- 
idence of overionization for Ar, given the respective ionization 
temperatures of kT z = 1.68 ± 0.07 and kT z = 1 .64+° ^. For Ca, 
the ionization temperature does seem higher than the contin- 
uum temperature: kT z = 2.0+° q 5 , with N H = 4.6 x 10 22 cirr 2 , 
and kT z = 2.04+°' 8 , with /V H '= 5.23 x 10 22 cirr 2 . However, 
we caution the reader that these estimates may not be reliable 
because the global spectrum originates from physically non- 
uniform regions of W49B that have different temperatures and 
abundances. A case in point in that the estimates of the elec- 
tron temperatures obtained from the global spectrum are signif- 
icantly lower those obtained from our spatially resolved spec- 
tral analysis of homogeneous regions (see Sect. 13.31 and ffii 
values in Table [3j. We therefore repeated our analysis for the 
more uniform central region indicated by a white ellipse in the 
lower left panel of Fig. [5] This is where the overionization ef- 
fects are expected to b e strongest in the scenario proposed by 
IKawasaki etafl d2005l) . Our analysis shows that these effects 
are not present. Using the same model described above in the 
2.75-6.0 keV band, we found for N H = 4.6 x 10 22 cirr 2 and 
Nh = 5.23 x 10 22 cm- 2 : respective electron temperatures of 
kT e = 2.4 + 0.1 keV and kT e = 2.25 ^f, Ar ionization tem- 
peratures of kT z = 1.9^2 an( ^ kT z = 1.8 + 0.1, and Ca ioniza- 
tion temperatures of kT z = 2.3_ ¥ ?:l and kT z - 2.2 + 0.2. Given 
that kT z is always slightly lower than, or consistent with kT e , 
we conclude that there is no indication of overionization in the 
central region of W49B. 

4. Discussion 

On the basis of low angular resolution ASCA line images, 
iFuiimoto et al.l ( 1 1 9951) argued for stratification of the SN ejecta 
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Table 3. Best-fit parameters for the spectra extracted in the seven regions shown in Fig.0](lower left panel), described with two 
thermal components in collisional ionization equilibrium. All errors are at the 90% confidence level. 

Parameter Region 1 Region2 

N a (l0 22 cm- 2 ) 4.5 ±0.1 4.22+°,"* 

kT x (keV) 0.90^-™ 0.75+°, : ° 3 

EM\ (10 20 cm- 5 ) 0.5+° 2 

kT 2 (keV) 3.0^ 2.19+° : ' 



EM;(10 20 cm- 5 ) 0.5 ±0.1 1.1 ±0.1 
Si/Si 2.2 ±0 .4 3.4+ 02 



J -0.2 ^- lv -0.06 

in-'" ,,,»-' i n 

!q = u.-t J -f_o,4 

S/5 4.1 ±0.5 4.1 ±0.3 

Ar/Ar Q 5.6% 3.8 ± 0.6 

Ca/Ca Q 8 ± 1 6.6+°! 

Fe/Fe 5.2+°, <j 5.1 ±0.3 

„2 
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4.5 ±0.1 
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4.61 ±0.05 


4 36+ 08 
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O Q9+0.06 
"■^^-0.08 


0.8 ±0.1 


0.96+°$ 


0.78+°$ 




1.2 ±0.2 


1 3+ 03 


0.8+° 3 




92+ ai 

"•^ -0.05 


2.6 ± 0.2 


9 oc+0.12 
' i --' J -0.09 


9 9+O.2 
'• -0.2 


i 77+O.03 
L - 1 '-0.09 


2.00 ± 0.04 




1 95+° 1 

7 -0.08 


1.4 ±0.3 


9 9+0.3 

•'••'-o.i 


1.3 ±0.1 


3.7 ±0.4 


2.9 ± 0.3 


9 4+0.3 
-0.5 


2.3 ±0.1 


2 r°' 3 


4 6+ 03 


3 6+ 03 


9 7+O.3 
^■'-0.5 


9 4+0.2 
'•^-o.i 


2.8 ±0.3 


a q+o.8 


4 1+0-3 


2.6 ± 0.6 


2.8 ±0.3 


3 4+ 03 


7 ± 1 


6 9+° 7 


4 7+0.8 
^■'-0.9 


4 4+0.4 
-0.5 


5 2+°- 5 

J " -0.7 


5 l+°' 4 


5 6+ 04 


2 6+ 02 

'■"-0.3 


2 5+ ai 

'• J -0.2 


1 6+ 02 



X z /d. o. f. 298.8/244 506.7/443 800.0/646 793.5/674 403.2/352 1059.3/844 566.3/434 



/l 




0.11*Jg 


0.08 ± 0.04 


0.08 ± 0.03 


10+ uue 

u - lu -0.05 


09 + 06 

"■"-'-0.04 


11 + 004 
"■ -0.03 


«i (cm 3 ) 




7!? 


7+4 

-2 


8±2 


5+ 3 

J -2 


7+4 

'-2 


5 5+ 3 


«2 (cirr 3 ) 


1-8^ 


"33 


2.6 ±0.2 


2 6+° 3 


2.3 ±0.2 


3 0+°- 8 

-'•"-0.2 


2.3 ±0.2 



* Emission measure per unit area. 



in W49B, with iron confined to the inner region and lower-Z 
elements (Si and S) forming an outer envelope. This hypothe- 
sis was supported by the analysis of the global spectrum, with 
the intensity ratios of Ly g and Heg giving d ifferent ionization 
ages for Si, S, Ar and Ca. lHwang et alJ ( EoO O) pointed out that 
while this conclusion is compatible with the AS CA spectrum, 
no difference in the ionization ages is required if the data are 
modelled with a two-component thermal plasma. The analy- 
sis of XMM-Newton data allowed us to discriminate between 
these two scenarios. Our spatially resolved spectral analysis 
has shown that a correct fit of the spectra requires two ther- 
mal components and that these two components are compat- 
ible with collisional ionization equilibrium, both at the center 
and in the outer parts of the remnant. Moreover, as explained in 
Sect. 13.21 our line images and equivalent width maps exclude 
a radial stratification of the elements in the ejecta, with Fe K 
emission confined to the inner part of the X-ray emitting re- 
gion and lower-Z emission lines originating in an outer shell. 
In fact, both the central barrel-shaped structure and the eastern 
elongated structure of the SNR have high abundances of Ar, Ca 
and Fe, while the western part has lower abundances. 

The physical conditions of the plasma are not uniform in 
W49B. The highest temperatures are found in the northern end 
of the elongated eastern structure (labelled a) and in the central 
barrel-shaped structure (labelled b). The lowest temperatures 
are found in the western part of W49B and in the south-eastern 
zone. Assuming pressure equilibrium between the two compo- 
nents, no significant variations in the plasma density across the 
remnant are observed. The observed X-ray emission originates 
in the SN ejecta, and Si, S, Ar, Ca, Fe and Ni are all highly over- 
abundant with respect to their solar values. They have almost 
the same values in structures a and b, and are lower (but still 
overabundant) in structure c. No clear evidence for X-ray emis- 
sion from the shocked ambient medium is present. In fact, the 
results of the spatially resolved spectral analysis suggest that 
the both the X-ray emitting components are associated with 
overabundant plasma. This is in agreement with the similar- 



ity of the morphology between the Si He-like emission blend 
(associated with the cooler component) and the other emission 
blends (associated with the hotter component). As explained 
in Sect. 13.31 however, the high absorption column inhibits the 
possibility of a detailed study of the soft component in that it 
is not easy to determine its morphology. Furthermore, the con- 
tribution of the hot component is significant even below 1.7 
keV, so that it is not possible to produce a clean image of the 
soft component. We therefore cannot rule out the possibility 
that the soft component may be associated with the shocked 
circumstellar or interstellar medium. 

It is interesting to consider to what extent the complex mor- 
phology of W49B, and the inhomogeneities in temperature and 
abundance, are caused by the dynamics of the SN explosion or 
by the remnant's interaction with its environment. Recent 2.12 
micron infrared imaging observations have revealed a shocked 
molecular H2 cloud that confines the remnant (Fig. [2}. The cir- 
cular distribution of the molecular cloud around the SNR has 
likely been shaped by strong stellar winds from the progenitor 
star pushing out the ambient interstellar material to its current 
location and forming a cavity. Infrared [Fe II] observations at 
1 .64 ymi reveal a barrel-shaped structure with four coaxial rings 
of warm gas, whic h has been i nterpreted as residuals of a strong 
stellar wind ( Keoha ne et all2 004). 

The morphology of the eastern region of W49B has clearly 
been shaped by the ambient molecular cloud. Despite their sim- 
ilar abundance values, the eastern a and central b regions have 
very different morphologies. The fact that region a is almost 
perpendicular to region b is probably due to its confinement by 
the large H2 cloud to the east. The [Fe II] observations show 
that the cavity is an inhomogeneous environment for the super- 
nova remnant expansion. 

In contrast, there are no observed interstellar structures that 
might have shaped the elongated structure of the central X-ray 
barrel b. This centrally barrel-shaped structure thus may have 
been produced directly by the supernova explosion or by pre- 
existing structures in the ambient stellar wind material. 
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The central barrel morphology seen in the Chandra Fe K 
image has been attributed to a bipolar jet, and taken as evi- 
dence that W49B is a y-ray burst (GRB) remnant 1 . The pres- 
ence of such an elongated morphology, however, does not 
necessarily reflect a jet-like explosion. A similar elongated 
central bright structure has be en observed with Chandra in 
the G292.08+1.8 remnant, but IPark et aD _2004_ showed that 
its emission is associated with normal chemical composition 
plasma that originates in shocked dense circumstellar material 
rather than shocked ejecta. Unlike for G292.08+1.8 SNR, our 
analysis suggests that the X-ray emission in the central X-ray 
barrel of W49B is dominated by ejecta, although we cannot rule 
out the possibility that the soft component comes from solar 
abundance material. The scenario proposed for G292.08+1.8 
can nevertheless be invoked for W49B. It is possible that the 
barrel-like X-ray emission reflects the morphology of the am- 
bient medium across the remnant. While no dense structures 
are visible toward the North and South of W49B (as shown 
in the infrared image of Fig. |3, we may assume the ambient 
medium to be enhanced in the central zone, forming a torus- 
like structure all across W49B. According to this scenario, the 
cold component may be associated with the thin belt of material 
heated by the transmitted shock travelling into the ISM struc- 
ture. This would explain the low values of the filling factor of 
the soft component. The hot component is associated with the 
ejecta interacting with the reflected shock. In this case, only the 
material expelled in the plane of the torus is interacting with 
the the reverse shock and heated by it. Notice, however, that 
according to our estimates of the plasma density for the spec- 
tral regions, the total mass of the ejecta in the hot component 
in these regions is ~ 8 M . This indicates that if we assume 
that we are observing only a small fraction of the ejecta we 
probably overestimate the total expelled mass. 

Another possibility is that the barrel is associated with a 
bipolar ejecta jet interacting with the cavity wall visible at the 
eastern end of the remnant. In this case we have to assume 
that the density contrast has driven a very strong reverse shock 
which has now reached the center of W49B . Moreover, we have 
found that structures a and b have similar abundance values and 
therefore, that a could be considered as the head of the eastern 
jet. The scenario that emerges is that this jet has been distorted 
to deviate mostly southwards by its impact with the H2 wall. In 
this framework, we would have a highly anisotropic bipolar jet, 
where the eastern arm is significantly more extended, hot and 
Fe-rich than the western one. 

Jet-like explosions have been investigated as possible 
mechanisms for y-ray bursts (GRB) because a spherically 
symmetric GRB would require an extremely large energy (sev- 
eral times the solar rest mass). A lin k between very ener - 
getic SNe and GRBs has been found bv lGalama et alJ (^998), 
who showed that the bright type Ic SN 1998bw (with an es- 
timated energy E ~ 3 x 10 52 erg. llwamoto et al Jl 19981) is the 
probable optical counterpart of the y-ray burst GRB 980425. 
Observations of SN 1998bw have shown st rong indications of 
asphericity in the explosion (Mazzali et al. 2001). Indications 
of aspherical explosions have also been found in many other 

1 http://chandra.harvard.edu/press/04_releases/press_060204.html 



hypernovae, i. e. supernova explosions with E <: 10 erg (for 
a review on hypernovae see lNomoto et all2003l) . 

The relationship between asymmetric hypernovae 
and GRB is in agreement with the collapsar model 
(MacF adven& Wooslevl 1 1999b . in which the collapse of 
a massive rotating core can generate a black hole and su- 
personic jets which propagate through the star. The nebular 
spectrum of SN 1998bw is well reproduced by such a model. 
iKhokhlov et alJ 0999) have shown that this mechanism 
produces a highly aspherical explosion with bipolar fast jets. A 
bipolar explosion does not necessarily imply a GRB, however, 
because a high explosion energy (^ 10 52 erg) is also required. 

To investigate the possibility that W49B is a GRB remnant, 
we compared the abundance values found in the a barrel (see 
Table 0} with those obtained from models of explosive nucle - 
osynthesis in bipolar explosions by Maeda & Nomotol (|2003). 
In particular, we compared the expected and observed ratios 
(X/X Q )/(Fe/Fe Q ) for X = S i, S , Ar, Ca, Cr, Mn and Ni (for 
Cr and Mn we assumed the same ratios as observed for Ar , 
in agreement with the estimates given in iHwang et alJ EoOO). 
None of the models can reproduce the entire observed abun- 
dance pattern. The hypernova models with zero age main se- 
quence progenitor mass Mzams = 40M o and explosion en- 
ergy E > 10 52 erg are strongly in disagreement with the ob- 
served values (with very large discrepancies for the Si, S and 
Ar abundances). We obtained better results with the less ener- 
getic models 25A (E = 6.7 x 10 51 and M ZAM s = 25M ) and 
25B (E = 0.6 x 10 51 and Mzams = 25M ). We also found 
relatively good agreement with the spherical explosion model 
25Sa (E = 1 x 10 51 and Mzams = 25M Q ). The predicted val- 
ues of these three models and our observational findings are 
summarized in Fig.[S] There is good overall agreement for Ar 
through Ni, but the observed abundances of Si and S are not 
well reproduced. 

While the observed abundance pattern in W49B does not 
favor an hypernova explosion of sufficient energy to generate a 
GRB, we cannot definitely exclude that W49B was related to 
such an event. 

5. Summary and conclusions 

Our study of the XMM-Newton observations of the Galactic 
supernova remnant W49B has shown that: 

- The remnant presents a complex X-ray morphology charac- 
terized by centrally barrel-shaped emission (b), terminated 
on the eastern side by an elongated perpendicular struc- 
ture (a) and on the western side by a more diffuse, nearly 
aligned structure (c). 

- The physical conditions of the plasma are not homogeneous 
throughout the remnant. Lower temperatures are found in 
remnant's western end, while hotter regions are located in 
the center and northeast. 

- A high Ni overabundance (Ni/Ni Q = 10+ 2 ) has been mea- 
sured in the central barrel, and the presence of Cr and Mn 
lines in the global spectrum is confirmed. 

- We found no evidence for overionization of the plasma in 
the center of W49B. 
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- The X-ray emission from all the observed regions of W49B 
originates in the ejecta, based on significant element over- 
abundances of Si, S, Ar, Ca and Fe. The equivalent width 
maps of W49B and a detailed spectral analysis of individ- 
ual regions both show that the central (b) and eastern (a) 
structures have similar abundances, and the western (c) re- 
gion exhibits lower abundances. This decrease is particu- 
larly pronounced for iron. 

- The X-ray morphology of W49B may be attributed to pre- 
existing structures in the ambient medium which may have 
produced an aspherical reverse shock, or alternatively to an 
aspherical supernova explosion in a wind cavity (bounded 
by molecular clouds on the eastern side of the remnant). 
Because it has abundances similar to those in the central 
barrel, the eastern structure a can be considered as the head 
of the jet. It is distorted and deviates mostly southward due 
to its impact with the dense molecular wall that confines its 
expansion to the east. 

- We investigated the possibility that W49B is a y-ray burst 
remnant by comparing the observed abundances with yields 
derived from supernova and hypernova nucleosynthesis 
models. We found better agreement with the less energetic 
(E ~ 10 51 erg) models. The association of W49B with a 
GRB is therefore not necessarily implied by the spectral 
data. 
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Figure 1. Left panel: Vignetting corrected EPIC count rate im- 
age of W49B in the 1-9 keV band. The image is adaptively 
smoothed (with signal-to-noise ratio 10) and background- 
subtracted. The 2 pc scale has been obtained assuming a dis- 
tance of 8 kpc. The count rate ranges between and 7.2 x 10~ 3 
s -1 . Right panel: Radio ima ge of W49B at 32 7 MHz (angular 
resolution 6"), obtained bv lLacev et alJ 1 1200 11) . We have super- 
imposed in black the X-ray contour levels derived from the 1-9 
keV EPIC image at 25%, 50% and 75% of the maximum. The 
color bar has a linear scale. 



Figure 2. Color composite Chandra X-ray image (in 
blue), infrared 2.12 jjm molecular hydrogen (in red), 
and 1.64 fim [Fe II] (in green) images of W49B (from 
http://chandra.harvard.edu/press/04_releases/press_060204.html) 
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Figure 3. pn (upper) and MOS (lower) global spectra of W49B 
in the 1-9 keV energy band, extracted from a circular re- 
gion covering the entire remnant (of radius 3.3' corresponding 
to ~ 7.7 pc). Line energies have been determined using nar- 
row gaussian components on top of a thermal bremsstrahlung 
continuum. Note that to reproduce the continuum at energies 
below the Si K lines, an additional bremsstrahlung component 
is required. 

Figure 4. Vignetting corrected and continuum-subtracted EPIC 
count rate images of the continuum (4.4 - 6.2 keV, upper left); 
S XV K (2.35 - 2.53 keV, upper right); Ar XVII K (3.05 - 3.2 
keV, central left); Ca XIX K (3.8 - 4.0 keV, central right); Fe 
XXV K (6.45 - 6.9 keV, lower left) and Fe XXV+Ni XXVII 
K (7.65 - 8.0 keV, lower right). All images are adaptively 
smoothed with signal-to-noise ratio 10. The 7 regions selected 
for spatially resolved spectral analysis are shown in the Fe 
XXV K line image, while the ellipse in the Fe XXV+Ni XXVII 
K line image indicates the region used for the detection of 
the Ni line (Sect. 13.4V The count rates range between and: 
8.6 x 10~ 4 s- 1 (continuum); 3.4 x 10~ 4 s" 1 (S); 1.7 x 10~ 4 s" 1 
(Ar); 1.9 x 10~ 4 s" 1 (Ca); 7.0 x 10~ 4 s" 1 (Fe); 2.1 x 10~ 5 s" 1 
(Fe+Ni). The color bar has a linear scale. 



Figure 8. Elemental abundances in the central X-ray barrel of 
W49B (red crosses, see Table 0}, normal ized to the Fe abun - 
dance relative to the solar values jAnders & Grevessd ll989). 
compared with the yields obtained from the aspherical explo- 
sion models 25a (boxes) and 25b (diamonds) and from the 
sphe rical model 25Sa (st ars) of explosive nucleosynthesis by 
Maed a & Nomotol |2003). The Cr and Mn abundances in the 
central barrel have been d erived from the Ar abundance fol- 
lowing lHwang et all (2000). 



Figure 5. Equivalent width maps of the S He-like and H-like 
line emission in the 2.35 - 2.7 keV band (upper left); of the 
Ar XVII K blend (3.05 - 3.2 keV, upper right); of Ca XIX K 
(2.35 - 2.7 keV, central left); of Fe XXV K (6.45 - 6.9 keV, 
central right) and of the Fe XXV+Ni XXVII K emission lines 
in the 7.65 - 8.0 keV band (lower left). In the lower right panel 
we have overlaid on the Fe XXV K equivalent width map the 
75%, 50% and 25% contour levels of the corresponding line 
image (Fig.|4] lower left panel). The ranges of the equivalent 
width maps are: 0.245 - 0.6 keV (S); 0.07 - 0.13 keV (Ar); 
0.085 - 0.18 keV (Ca); - 12.1 keV (Fe); 0.315 - 1.4 keV 
(Fe+Ni). 



Figure 6. Vignetting corrected and continuum-subtracted EPIC 
count rate image of the Si XIII Ka (1.77 - 1.9 keV, left panel) 
and Si XIV Ka (1.96 - 2.06 keV, right panel). The images are 
adaptively smoothed with signal-to-noise ratio 10. The count 
rates range between and 2.3 x 10~ 4 s _I for Si XIII and 1.7 x 
10~ 4 s _1 for Si XIV. The color bar has a linear scale. 



Figure 7. pn (upper) and MOS (lower) spectra of region 2 with 
the corresponding best-fit model and residuals obtained with 
one thermal component (left panel) and with two thermal com- 
ponents (right panel). The contribution of the two components 
for the MOS camera is indicated in the right panel in cyan for 
the cooler component and blue for the hotter component. 



